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O Abstract 

■ Solar neutrino physics enters a stage of precision measurements. In this connection we 

present a precise analytic description of the neutrino conversion in the context of LMA 
MSW solution of the solar neutrino problem. Using the adiabatic perturbation theory 
we derive an analytic formula for the v e survival probability which takes into account the 



non-adiabatic corrections and the regeneration effect inside the Earth. The probability is 
, averaged over the neutrino production region. We find that the non-adiabatic corrections 

J> ' are of the order 10~ 9 — 10~ 7 . Using the formula for the Earth regeneration effect we 

discuss features of the zenith angle dependence of the v e flux. In particular, we show that 
, effects of small structures at the surface of the Earth can be important. 

o 

^ ; 1 Introduction 



The LMA MSW solution [1, 2] has been identified [3] — [15] as the correct solution of the solar 
neutrino problem. The 2z/ conversion probability of this solution gives a very good description 



I 

D \ of all available data: no statistically significant deviation has been found so far. New physics 
^ I effects beyond LMA, if exist, are below few per cent. 

The program of future solar neutrino studies includes 
^ ■ 1). further tests of the LMA solution, in particular, searches for signatures of this solution 
such as the Day- Night asymmetry and the distortion ("upturn") of the boron v e spectrum at 
low energies; 

2) . precise determination of the oscillation parameters, especially the 1-2 mixing angle; 

3) . searches for the sub- leading effects which originate from 

- 1-3 mixing, 

- sterile neutrino mixing, 

- non-standard neutrino interactions, 

- spin-flavor flip in the magnetic fields of the Sun, 

- violation of the fundamental symmetries (CPT, equivalence principle, etc.). 

Already the present solar neutrino measurements have sensitivity at the level of few per 
cent. For instance, the predicted day- night asymmetry of the SuperKamiokande signal is about 
2% which is comparable with the existing la experimental error [6]. At SNO one expects the 
2 — 4% asymmetry, consistent with the experimental result [3] at the lcr level. 

Future experiments will have substantially higher sensitivity [16, 17, 18, 19]. The solar 
neutrino studies enter a phase of precision measurements. 
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In this connection it is important 



• to give precise description of the LMA conversion, both in the Sun and in the Earth, 
taking into account various corrections; 

• to estimate accuracy of the approximations made; 

• to find the precise analytic expressions for probabilities and observables as functions of 
the oscillation parameters (Am 2 , sin 2 612) ■ This will help to test the LMA solution and 
to search for physics beyond LMA. 

We address these issues in the present paper. In section 2 we consider the non-adiabatic 
corrections to the LMA conversion probability. We calculate these corrections for propagation 
inside the Sun and the Earth. In section 3 we obtain the analytical formula for the probability 
averaged over the distribution of neutrino sources. In section 4 we derive the analytic formula 
for the u e regeneration effect in the Earth. We present our conclusions in Section 5. In the 
appendices A and B, alternative derivations of formulas for the regeneration factor are given. 

2 Non-adiabatic corrections to the LMA solution 

According to the LMA MSW solution of the solar neutrino problem, a conversion of the solar 
electron neutrinos is driven by mixing of the two active neutrinos, ^/ = (u e , z/ a ) T : 

* f = u(e)* maaa , (i) 

where, in general, the mixing matrix is determined as 

U(a)=( C ° Sa 8ina V (2) 



sm a cos a 



and ^f m ass = {yii ^2) is the vector of mass states. 
2.1 LMA and Adiabaticity 

The main feature of the LMA solution is the adiabaticity of conversion. According to LMA the 
averaged 2v survival probability of the electron neutrinos is given by the adiabatic formula [2, 
20, 21]: 

P ee = ^(1 + cos 20^ cos 20). (3) 

Here 6 is the vacuum mixing angle, 6® n = 9 m (x ) is the mixing angle in matter in the neutrino 
production point, xo, and the mixing angle in matter is determined by 

Here Am 2 is the mass squared difference, E is the neutrino energy, V is the potential, Gp is 
the Fermi coupling constant and n e (x) is the number density of electrons in the point x. 
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How precise the expression (3) is, and what are the non-adiabatic corrections? To answer 
these questions, we will elaborate on the adiabatic perturbation theory. 

Dynamics of neutrino conversion is described in terms of the instantaneous eigenstates 
of the Hamiltonian in matter, v\ m , V2m- Representing an arbitrary neutrino state as \v) = 
ipimWim) + ^2m\ v im) i we can write the evolution equation in the base (vi m , z/ 2m ) as [2, 20, 22] 

.d_(4>im\ = ( -^nr -iO m (x)\ f^ lm \ 

where 



A(x) = Am 2 ^(cos 29 - 2EV(x)/Am 2 ) 2 + sin 2 29 (6) 

and 

a d6 m (x) _ EAm 2 sm29dV(x) 

m[x) = ~dx~ ~ A{xf dx~- { ) 

The adiabatic approximation corresponds to a situation when 

7-^«l, (8) 

and the off-diagonal terms in the Hamiltonian (5) can be neglected. In this case there are 
no transitions between the eigenstates, and therefore the eigenstates propagate independently. 
The solution of (5) is straightforward: 

^(x) = S ad (x,x )V m (x ), (9) 

where 

•*->-(£«). *.<*»> (10) 

and the adiabatic evolution matrix is 

S ad (x,x ) = S ad ^)=[ e Q J Hx) ). (11) 
The adiabatic phase $(x) equals 



<$>{x) = j E jjx'A{x'). (12) 



A state initially produced as the electron neutrino, ipim(x ) = cos#^, ip2m( x o) — sm $m> 
evolves as 

v{x) = cos9° m e^ x K lm + sme° m e-^u 2m . (13) 

The incoherent survival probability (3) can be immediately obtained by averaging the matrix 
element squared, \{u e \u{x))\ 2 , over the oscillations. 
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2.2 Non-adiabatic corrections 



The non-adiabatic corrections correspond to transitions between the instantaneous eigenstates. 
We calculate these corrections by solving the equation (5). We will implement a perturbation 
theory using the fact that for LMA the adiabaticity parameter (8) is very small. 
Let us search for the solution of the equation (5) in the form 



- c *(x) i ) \wi(x) 

e i*(x) c^e-^ 



ipim(xo) 

^2m(x ) 



(14) 



where \c(x)\ <C 1 is supposed to hold everywhere along the neutrino trajectory. (We check this 
a posteriori). 

The expression (14) can be rewritten as 



V m (x) = S(x,x )ty m (x ), 

where the evolution matrix equals 

S{x,x Q ) = CS ad ^ 



e i<s>{x) c{x)e~ i ^^ 
c*(x)e^ e - **^) 



(15) 



(16) 



Inserting (14) into (5) we find the differential equation for c(x) from the condition that the 
off-diagonal elements of the evolution equation for ip im are zero: 



l Tx c{x) 



A(x) 
^2E 



c(x) - iO m {x). 



(17) 



Here the first order terms in c(x) and 9 m are kept only. In this approximation the energy gap 
between the states coincides with the adiabatic split A(x) given in Eq. (6). 
The solution of equation (17) can be written in the following form 



c(x) 



ax' 



.-co 



-' / 

J X 



IE 



The integration constant is fixed by the condition : c(x) — > as 9 m — > 0, so that c(xq) = 0. 

Since for the LMA solution the phase <&(#) is a fast oscillating function, the integral in (18) 
can be calculated using the following formula (essentially the integration by parts) 



t f(x)e l9{x) dx 

J a 



J(x) , f(x) f(x)g"(x) 



g'yx) 



g' 2 (x) 



g' 3 (x) 



+ 0(l/g") 



(19) 



which is valid for smooth functions of f(x) and g(x). Here g'(x) = dg(x)/dx and f'(x) = 
df(x)/dx. The formula gives very good approximation if f(x)/g'(x) <C 1. In the case of inte- 
gral (18) this condition coincides with the adiabaticity condition (8) which is well satisfied. 
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According to (19) and (18) we find 



.IE d6 m (x) 

c(Xf) — —i — exp 

v 11 A(x) dx 



- f 

Jx f 



dx 



AW 

2E 



X=Xf 
X=XQ 



= -*sign(# m )^exp [-2i($(x) - 



(20) 



X=XQ 



or explicitly 



, , 2E 2 Am 2 sin29dV(x 

c \ x f) = ~ % a 7Z\a 7T^ ex P 



dx 



Jxf 2E 



X=XQ 



(21) 



where sign(0 m ) = 9 m /\9 m \. 

We will apply this formula for neutrino propagation inside the Sun in section 2.3 and inside 
the Earth in section 4. 



2.3 Non-adiabatic corrections for propagation inside the Sun 

The survival probability with the adiabaticity violation effect included can be written as 

1 + (1 -2P C ) cos 2^ cos 29 



(22) 



where P c = \c(xf)\ 2 is the jump probability — the probability of transition z/ 2m — > u lm on the 
way from x to Xf. 

Let us calculate P c . Notice that for the LMA solution, one can not use the Landau-Zener 
formula [23] for P c for the following reasons. The mixing angle in the final point of evolution, Xf, 
is large. The adiabaticity parameter \AE9 m /A\ is of the same order for all points inside 0.3 of the 
solar radius [24]. The point of maximal adiabaticity violation is not the resonant point, though 
not far from it. Moreover, the resonance layer defined as \2EV — Am 2 cos2#| < Am 2 sin 29, is 
broad since the vacuum mixing angle is large. Futhermore, significant part of the neutrino flux 
is produced inside the resonance region or does not cross the resonance region at all. 

The double exponential formula [25] is not valid too. It requires production of neutrinos far 
above the resonance region in the density scale. This formula is not applied in the range 



Am 2 
~2E 



cos20~(1.6-8.O)-lO" 



eV 2 
MeV 



(23) 



for which the density in the production point turns out to be close to the resonance density. For 
the best fit values of the LMA oscillation parameters the range (23) corresponds to E — (2 — 15) 
MeV, that is, to the region of interest. 

Let us apply the results of section 2.2 for calculation of the non-adiabatic corrections. Notice 
that at the surface of the Sun the effective potential V is negligible and 9 m can be taken zero. 
Using Eq. (21) we find the transition amplitude c(xf) on the way from the production point 
to the surface of the Sun in the leading order approximation as 



c(x f ) 



2E 2 Am 2 sin 29 dV(x) 



A(xf 



dx 



x exp 



X=XQ 



x f , A(x) 
i I dx-^-f 
x 2E 



(24) 



5 



K 


pp 


8 £? 


VA N 


15 


17 F 


7 5e 


pep 


hep 


V^(10- 12 eV) 


4.68 


6.81 


6.22 


6.69 


6.74 


6.16 


5.13 


3.96 




0.109 


0.010 


0.054 


0.013 


0.012 


0.029 


0.076 


0.165 



Table 1: TTie average value of potential Vk and the corresponding value of AV^/V^ for different 
component of the solar neutrino spectrum. 



Then the probability of non-adiabatic transition is given by 



Pc = \c(Xf)\ 2 = 



7 2 (zo) 



1 11 



16tt 2 h 2 (x ) 



2EV(x )Am 2 sin 26 



where 



h(x) = V(x) 



dVjx) 
dx 



losci^C) 



A(x ) 



AitE 
A{x) 



(25) 



(26) 



are the density height and the oscillation length in matter. 

The transition probability P c (25) depends only on parameters of the initial (production) 
point. One can understand this by noting that l OS c{ x ) ^ h c (x). Therefore many oscillation 
lengths are obtained on the distance in which the potential changes sizably. The non-adiabatic 
corrections are averaged out being negligible along the trajectory of the neutrino except for the 
boundaries of trajectory, i.e., around the production point or the point at the surface of the 
Sun. At the surface of the Sun the contribution can be neglected because the potential is zero. 

The probability is determined basically by square of the ratio of the oscillation length and the 
density height. Second factor in (25) is of the order one. So, essentially P c < [l osc (x ) / 47ih(x )] 2 . 
Using the best fit values of the LMA oscillation parameters we find from (25) 



Pc = (10 



-9 



10 



-7\ 



E v 



10 MeV 



) 



(27) 



Here the numerical prefactor depends on the production point. As a function of Xq, the prob- 
ability P c reaches maximum at around (0.1 — O.2)i? , where the potential doesn't drop down 
substantially and h(x ) reaches its almost minimal value due to increase of the density gradient. 
The corrections are negligible in the whole relevant range of neutrino energies and production 
points. The probability (27) strongly differs from what one would get using double-exponential 



formula 



^ e -4nh/l osc < 10" 400 . 



Notice that the jump probability (25) equals (up to factor 4) the adiabaticity parameter 
in the production point squared, as is expected in the adiabatic perturbation theory. This 
contrasts the Landau-Zener probability, P c ~ exp(— n/2j), which is essentially non-perturbative 
effect. 

3 Averaging over production region: analytic results 

In the adiabatic approximation the survival probability depends on the potential (density) in 
neutrino production point r: 



Pee = Pee(V ), V = V(x = r). 



(28) 
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Observables at the Earth are determined by the survival probability averaged over the 
neutrino production region: 

^ = JdrGKpP^r)_ = pp, pep, Be, N,0,F,B, hep, (29) 

fdr G K (r) 

where Gk(j') is the distribution of sources of the K component of neutrino spectrum. The 
distributions are different for different components. 

Let us introduce the average value of the potential in the production region for the type K 
neutrinos: 

___ fdr G K (r)V(r) 
VK = IdrG K (r) ■ (30) 

We will use the fact that in the effective production region, V(r) deviates weakly from Vk- 
Therefore the survival probability can be expanded in series around Vk' 



P ee = P ee (V K )+ (y-V K ) + -[—g) {V-V K f + ... (31) 




Inserting this expression into (29) and using the definition (30) we find 

P p (v \ 3E " sin 2 2gcos2g(cos2g-2^/Am 2 ) 

eeK K) (Am 2 ) 2 [(cos2# - 2EV K /Am 2 ) 2 + sin 2 29f/ 2 K ' 1 J 

where 

2 _ fdrG K (r)(V(r)-V K y 
Vk ~ fdr G K (r) • (33) 

Notice that the correction to P ee appears in the second order of the deviation of potential from 
the average value. The expression for probability (32) can be rewritten as 

P K = ± + ±(1-6 K ) cos29 m {V K ) cos 29, (34) 

where the correction 5k equals 

3 (2ffV^/Am 2 ) 2 sin 2 2fl AVg 

K 2[(cos2^-2E^/Am 2 ) 2 + sin 2 2^] 2 V K ' 1 ' 

In the Table 1 we present the average values of potentials and the corresponding second order 
deviations from the average values for eight types of solar neutrinos. We use the distributions 
of neutrino sources from the BP2000 model [26]. The expansion parameters AV K /V K are all 
small, especially for the boron neutrinos which have the narrowest distribution of sources. 

In Fig. 1 we compare the probability Pk obtained from the approximate analytic formula 
(34) with results of numerical calculations, P' K . Pk/ P'k — 1 is shown as a function of E j 'Am 2 . 
The lines have been cut for K = pp, Bel ', pep, N, O, F because of their lower energies in com- 
parison with the energies of hep and boron neutrinos. The plot shows that the analytic formula 
is rather precise. In particular, the deviations are extremely small ( < 10~ 3 ) for small and 
large values of E/Am 2 . Relatively large deviations can be seen in the intermediate region of 
E/ Am 2 . For example, for K = hep the magnitude of Pk/P'k ~ 1 reaches maximum 1.8% at 
around E/Am 2 m 34 x 10 10 eV _1 . The corrections 5k are important: e.g., for the hep neutrinos 
the deviation would be up to 6% without 5k ■ 
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0.01 




Figure 1: Deviations of the probability Pk given in formula (34) from the numerically calculated 
probability, P' K for different components of the solar neutrino spectrum. 

4 The Earth matter effect: analytic study 

The solar neutrinos arrive at the surface of the Earth as incoherent fluxes of the mass states. 
The mass states oscillate in the matter of the Earth producing partial regeneration of the 
electron neutrino flux [27, 28, 29, 30, 31, 32, 33, 34]. Previously the effect has been described 
in one or two layers approximation. In the later case, interference effects of contributions from 
the core and the mantle have been discussed [33, 34]. In this section we will study effects for 
the realistic density profile of the Earth. 

4.1 Regeneration factor and the Earth density profile 

The probability of — ► v e transition can be written as 

P{y 2 -> v e ) = sin 2 9 + f reg , (36) 

where f reg is the regeneration factor which describes the Earth matter effect. In the absence 
of matter (i.e., during the day) f reg = 0. Using the definition (36) we find the v e survival 
probability with the regeneration effect included as 

Pee = ~(1 + cos 26° m cos 26) - cos26° m f reg . (37) 
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r : radial distance 

x : trajectory coordinate 

R E : radius of the Earth 

R. j : outer radius of the 
ith shell (i=l,...,9) 
Vj : potential of the ith shell in 

R i <r<R i _ 1 (i=l,...,9,R 9 = 0; 



L =L : total length of trajectory 



in Earth 
L. : length of trajectory 

for r<R (i=l,...,8) 



Figure 2: Structure of the Earth density profile. We indicate notations used in the text. 



Notice that the mixing angle (P m in the neutrino production point in the Sun determines the 
mass {yi, z/ 2 ) composition of the neutrino flux which arrives at the Earth. 

The essential feature of the LMA solution is that the Earth matter effect is small. This 
smallness is characterized by the ratio 



7/ 



2EV 
Am 2 



(38) 



where is the Avogadro number. We will use 7] as the expansion parameter. 
The i>2 — > v e transition probability can be written as 

P(y 2 _> Ve ) = \(u e \U(9 mR )S(x f ,x )U\9 mR )U(9)\u 2 )\ 2 

= \(ise\u(e mR )s(x f ,x )tf(e mR -e)\v 2 )\ 2 , 



(39) 



where 8 mR is the mixing angle in matter at the surface of the Earth and the matrix S(xf,xo) 
given in (16) describes evolution of the neutrino eigenstates in matter. Noting that in matter 
of the Earth 

EVn 

sm{9 mR - 6) w — | sin 29 < 1, (40) 
Am 2 

we find from (14), (36) and (39) an expression for the regeneration factor in the lowest order 
in c(x) and sin(# m # — 6) as 



reg 



^| sin 2 29 sin 2 $(x f ) + sin 29Re{c(x f )}, 



(41) 
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where V R is the potential at the surface of the Earth and <&(£/) is the total phase acquired 
along the trajectory in the Earth. 

For the profile with slowly changing density (lowest adiabatic approximation), c ~ 0, we 
obtain 

The prefactor (the depth of oscillations) is determined by the potential at the surface of the 
Earth, whereas the phase is given by the integral along whole trajectory. For one layer with 
constant potential (density), and therefore c = 0, the regeneration factor (41) or (42) is reduced 
to the well known expression: 

Here L is the distance traveled by neutrino in the Earth and l m is the oscillation length in 
matter. 

Let us consider a neutrino propagation in realistic density profile of the Earth. The profile 
can be described by n nearly spherical shells of matter with sharp (step-like) density changes 
at the borders of shells and slow variation of density in layers between the borders. According 
to the PREM model n = 9 [35]. So, in ith shell (i = 1, • • • ,n), the potential Vi is a smooth 
function of the radial distance r. Crossing j shells corresponds to crossing 2j — 1 layers (see 
Fig. 2). We denote by Ri_i the outer radius of ith shell, so that R corresponds to the radius 
of the Earth: R = R E . 

The trajectory of the neutrino is characterized by the zenith angle 0z- We determine a 
position of neutrino along trajectory by the coordinate x with origin in the center of trajectory, 
so that 

x e [-L/2, L/2], x 2 = r 2 - R% sin 2 6 Z . (44) 

Here L is the total length of the trajectory in the Earth. The length of the part of trajectory 
inside border Ri is given by 

Li = yjR 2 -R 2 E sm 2 9 z , (45) 

and L = L by definition. 

We introduce the adiabatic phase <3>j acquired by neutrinos in the interval —Li/2 < x < L { /2 
along the trajectory, that is, inside the outer border of the (i + l)th shell: 



W2 A(x) 
ax - - 

U/2 AE 



j (Lr 

"Li/2 

dx 

Li/2 



/. 



Am 2 1 nmr . , Esm 2 26^ 2/ , 



0,---,n-l. (46) 



A(x) is given in Eq. (6). Here we keep the order V 2 term since due to integration its contribution 
to the phase is not negligible. 
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At the borders of shells there are jumps of the potential and hence the discontinuities of 
the mixing angle in matter. We denote them as 



AVi ^Vi+^Rj-ViiRi), i = 0,..-,n-l, (47) 

A9 mi =0 m (^i+iW)-MW)), i = 0,.--,n-l. (48) 
At the surface of the Earth, we obtain 

AV = V R , A9 m0 = 9 mR -9. (49) 

Furthermore, we find 

smA9 mi w ^ — ^- sin 20, cosA0 mi f« 1, i = 0, • • • , n — 1. (50) 
Am 2 

Corrections to (50) are of the order if, and hence negligible. 

Smooth variation of the potential in each shell of the Earth can be approximated by the 
analytic formula [36]: 



Vi = V A 



2 4 



p2 1 " p/ 



n. (51) 



Let us find an analytic expression for the regeneration factor using the density profile de- 
scribed above. The problem can be solved in two steps: (1) computation of the non-adiabatic 
corrections to propagation within a given layer, Aff, (2) computation of effect of the borders 
between layers, A// ump . So that 

(52) 

i i 

The virtue of the LMA solution is that it enables us to study both effects using the same 
formalism of the adiabatic perturbation theory. 



4.2 Non-adiabatic corrections in a layer of the Earth 



Let us compute the non-adiabatic corrections for one layer. Suppose a neutrino trajectory 
crosses the iih layer with the borders at x — Li/2 and x = Lj_i/2. Using (44) and (51), the 
potential in this layer can be expressed in terms of the trajectory coordinate as 



Vi = V A 



x 



x 

R% 



Here 



a. 



ai + fa sin 2 9 z + 'y i sm 4 9 z , f3- = fa + 2% sin 2 Z , l[ = li- 



According to (45) and (54) the gradients of potential at the borders equal 



dVAx) 



dx 

dVjjx) 
dx 



-hi. 

" 2 



2V A 
Re 

2Va 
Re 



y 2 -sin 2 9 z (fa + 2 liy 2 ), 



y 2 _ 1 -sm 2 9 z (fa + 2^yl 1 ). 



(53) 

(54) 

(55) 
(56) 
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where i/i = Ri/R E - Then for this layer Eq. (21) gives the amplitude of non-adiabatic transition 



= -lie 



:2m E 2 Am 2 sm29dV(x) _ 2m 



x) 



A 3 (x) 



AiE 2 sin 29 
(Am 2 ) 2 R E 



V A e^ 



dx 



{(3 i + 2 li y 2 _ 1 )Jy 2 _ l -sm 2 9 z e 



((3 l + 2^y 2 )Jy 2 -sm 2 e z e-^ 



(57) 



where phases $j are defined in (46). Inserting this expression into (41) we obtain the non- 
adiabatic correction from this layer, A/j to the regeneration factor as 

4£ 2 sin 2 2fl 
A/i = — — , XOn Va 



(58) 



^ , , . t (A + 2Ti3/,-_i) V3/i-i - sin 2 # z sin($ - 
yAm^y R E L 

-(A + 2 ll y 2 )^y 2 -sm 2 ^sin($„ - $j) 
The ratio of the absolute value of correction and the adiabatic term equals 



i f ad i 

I J reg I 



2E 



V A ((3i + 2 llV 2 



Am 2 sm29R E 



V, 



R 



y 2 - sin 2 Z ~ 



2E 1 



(59) 



where i?s plays the role of typical scale of the density change. As an example, let us consider the 
layer between 0.895^ and 0.5A6R E . In this layer a = 3.156, (3 = -1.459 and 7 = 0.280 [36]. 
From (59) we obtain that the non-adiabatic correction to the regeneration factor is about 
(1 - 2)% at E = 10 MeV. 

Notice that for some particular values of energies and 9 Z , the contributions A/j from different 
layers % may sum up "constructively" producing larger effect. In this connection let us notice 
the following. 

1) The enhancement effect may occur for exceptional values of E and $z and therefore any 
realistic averaging over E and integration over 6 Z will wash it out. 

2) The enhancement can not be large (proportional to the number of layers, n) since (i) 
only few layers give significant contribution and for the rest, the effect is below 1%; (ii) there 
is a systematic cancellation of contributions from the upper limit of integration in A/j and the 
lower limit of integration in A/j_i (the adiabatic phases are the same for both); (iii) typically, 
contributions from two layers of the same shell have opposite signs. 

So, we conclude that the non-adiabaticity within layers of the Earth can be safely neglected. 



4.3 Effects of several layers 

The jumps of potential between the layers strongly violate the adiabaticity and on the first 
glance, the adiabatic perturbation theory can not be applied. We show, however, that the 
results for non-adiabatic case obtained in section 2.2 can be also used here. The key point is 
that for the LMA parameters the Earth matter effects are small, whatever the density profile in 
the Earth is. As a consequence, variations of the mixing angle in matter are small: |A0 m | <^ 9, 
and essentially the expansion parameter here is rj. 
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Figure 3: The regeneration factor as function of cos 6 Z for E = 10 MeV, Am 2 = 6.3 x 10" 5 eV 2 , and 
tan 2 8 = 0.4. We compare result of numerical computations with analytic result (63) for the PREM 
model. 



Consider a neutrino trajectory which crosses In — 1 layers (n shells). In the points of the 
trajectory x = =F-Lj/2 neutrinos cross the borders of shell with r = Ri, as illustrated in Fig. 2. 
The corresponding potential jumps equal ±AT^ for x = =pLj/2. Using the expression (7) we 
obtain in the lowest approximation 



e m (x) = ^-^ E ^ 



i=l 



5[x + 



L, 



Six 



Li 



(60) 



As it has been shown in section 4.2 to a good approximation one can take 9 rn = everywhere 
outside the borders. 

The evolution equation (5) can be averaged in small intervals Ax 1/A(x) to eliminate 
5-functions which originate from 9 m . However, this is not necessary since in the expression for 
c(x) in (17) 9 m is integrated anyway. 

Plugging expression (60) into Eq. (18) we obtain 



c(L/2) 

where $j are defined in (46). 



Esin29 
Am 2 



n-l 



(61) 
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Inserting the real part of c(L/2) into (41) we find the regeneration factor for the case of n 
shells crossing: 



freg = sin 2 26 sin 2 $ + J2 sin2 26 sin $ * sin ( 62 ) 

In the Appendix A we present the rigorous derivation of this factor considering evolution in 
the sequent layers explicitly. The results of two approaches coincide exactly in the first order 
in £17 Am 2 . 

Noting that Vr = AVq is the jump of potential at the surface of the Earth, we can rewrite 
the expression (62) in the following compact form: 

2E sin 2 29 n ~ 1 

freg = A ^ 2 sin $ £j AV; sin $i. (63) 

So, f reg is proportional to the sum of similar terms which correspond to the borders of the 
shells. Each term is the product of the potential jump at a given border and sine of the total 
adiabatic phase acquired on the part of trajectory inside a given border (that is, from —Li/2 to 
Li/ 2 for the border i). The sum runs over all borders including the surface of the Earth. The 
expression (63) corresponds to the symmetric density profile. The zenith angle dependence of 
the regeneration factor appears via the phases: $j = &i(0z)- 

The formula (63) (which is the main result of our study) allows us to get complete under- 
standing of the Earth matter effects including effects of complicated shell structure. Apparently, 
this is not possible using the one layer approximation (42), where the interference terms induced 
by different shells are absent. 

In Fig. 3 we compare the zenith angle dependence of the regeneration factor computed 
using the analytic formula (62) with the one obtained by the exact numerical integration for 
the PREM profile. Two results coincide extremely well. One can see that the analytic formula 
reproduces quite precisely the magnitude and the phase structure of the regeneration factor. 

Let us mark some features. The change of the oscillatory behaviors for cos 9z > 0.83 is 
induced by the sharp density jump at the border between the mantle and the core of the Earth, 
at r = 0.5AR E - Notice that at cos6>^ > 0.83 the amplitude of oscillations for some periods 
increases, however, the frequency of large peaks becomes lower. So that the average value of 
f reg does not increase in comparison with the value for cosdz < 0.83. 

For small cosdz the dependence of f reg on 9z is a result of interference of terms in (63) which 
correspond to the outer shells of the Earth. To understand this dependence, it is convenient to 
introduce the phase ipf. 

<Pi = $ - (64) 

so that ipi/2 is the phase acquired by neutrino on the way from the surface of the Earth to 
r = Ri. Using we can rewrite the expression for regeneration factor (63) as 



2E sin 2 29 2=J _ r » , '< 



Jreg A 2 1^^ V * 

lA " b i=0 



sin 2 $ cos ¥i sm 2^o si n V 9 ? 



(65) 



Apparently, if averaging over <fi occurs only the term with % — survives (ipo = 0) in f reg which 
is reduced to the adiabatic expression for one layer. 
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Figure 4: 

E = (8- 



The regeneration factor averaged over the energy intervals (a) E = (9.5 — 10.5) MeV; (b) 
11) MeV. For oscillation parameters we take Am 2 = 6.3 x 10~ 5 eV 2 and tan 2 6 = 0.4. 



The increase of regeneration factor with cos9z in the range 0.2 — 0.5 is related to the effect 
of three density jumps near the surface of the Earth: according to the PREM profile they are 
situated at depths 10 km, 22 km and 31 km correspondingly. The distance, L — L iy on which 
<fi is acquired depends on the zenith angle as 

Re — Ri 



L-Li 



cos 0; 



cos 9 Z > Jl - Rf I R% 



(66) 



(Re — Ri is the depth from the surface of the Earth to the borders of the ith shell). In the 
case of small cos8z (cos8z ^ 0.2), the distance L — Li can be of several hundreds kilometers 
which is comparable to or larger than the oscillation length. Furthermore, L — Li and ipj are 
fast changing functions of 9 Z . So, are large and different for different i. Therefore, the terms 
sin $ sin $j for different % (i = 1,2,3) are quite different and therefore partially cancel each 
other ("interfere destructively"). 

In constrast, for cos#z > 0.5, the distances L — Lj (i = 1, 2, 3) for the outer shells become 
much smaller than the oscillation length and they slowly change with cos#z. In this case the 
phases (fi (i = 1, 2, 3) are all small and sin $ sin <3>j sin 2 $ . So, for cos 9 Z > 0.5 the effects of 
outer shells "interfere constructively" producing larger regeneration factor. It is then possible 
to account these close layers effectively as a single layer, as it was done in Ref . [36] . Increase of 
the regeneration factor in the transition region cos#z = (0.2 — 0.5) corresponds to convergence 
of the term sin $o sin <E»j (i — 1, 2, 3) to sin 2 <3> - 

4.4 Averaging over the neutrino energy 

For the LMA solution the oscillation length in the Earth is small: l m ~ l v Re- Since the 
time of the neutrino detection is well known, averaging over the zenith angle can be avoided, 



15 



and in fact, in the unbinned analysis of the data developed recently [6] one needs to know the 
zenith angle dependence without averaging. At the same time since the recoil electron (and 
not neutrino) energy is measured and a detector has finite energy resolution, averaging over 
the neutrino energy occurs. 

In the leading approximation the phase equals (pi Am 2 (L — Li)/(AE). Therefore the 
energy resolution AE corresponds to averaging over the interval of phases: 

AE 

A(fii w <Pi-jjr- (67) 

If L — Li ^> l m , so that ipi >> 1, the interval Aifi can be large, thus leading to strong averaging 
of terms sin <3> sin <3>j in (63). This happens to the contributions from structures situated far 
from the surface of the Earth. In Fig. 4 we show the result of averaging of the regeneration 
factor folded with the cross section of the neutrino-electron elastic scattering over two different 
energy intervals. Comparing with Fig. 3, one sees that the complicated oscillatory pattern 
produced by the density jumps in the central regions of the Earth is strongly averaged when 
cos 9 Z is large (for general analysis of this effect see [37]). According to Fig. 4 for cos6 z > 0.4 
the regeneration factor oscillates with small depth around f reg ss 1.5%. This happens because 
the main term sin 2 $o is strongly averaged too. 

In contrast, for cos 6z ^ 0.2 only the outer structures of the Earth can contribute and 
the averaging is not as efficient as for large cos9z- Indeed, for the borders of outer shells 
(i = 1, 2, 3) the distance where Lpi is acquired, L — Li < 2Re\Ji — Rf / R%, can be about several 
hundreds kilometers, that is, comparable with the oscillation length in matter. In this case 
the phase interval A<pi ps (piAE/E is still not large enough to give sufficient averaging. For 
cosOz < 0.4, A$o ~ &oAE/E is also small and averaging is weak. Furthermore, in the interval 
cos#z = (0.2 — 0.5) the regeneration factor increases with cos8z- The reason is that in this 
interval the main term, sin 2 $o, starts to "interfere constructively" with the terms produced by 
the outer shells, sin $ sin $j (i = 1, 2, 3), as it has been discussed in section 4.3. 

In Fig. 5 and 6 we show dependence of the regeneration factor for the charged current 
events at SNO. Here we have taken into account the energy resolution of the SNO detector 
and also performed integration over various energy bins of the observed kinetic energy. On the 
basis of our analytic formulas and discussion, the interpretation of results of Fig. 5 and 6 is 
straightforward. 

4.5 Small scale structures: general density profile 

There are small scale structures in the outer mantle of the Earth of depth (~ 10) km in which 
matter has quite different densities (e.g., ocean, rock and soil). In contrast to the ideal PREM 
model, these structures are not isotropically distributed and can be quite complicated. 

In section 4.3 we have shown for the ideal PREM model that contributions produced by 
structures close the surface of the Earth interfere destructively for small cos 8z- Furthermore, 
averaging over the energy doesn't smooth the dependence of these contributions on cos Qz com- 
pletely. This produces an uncertainty for the future high statistics solar neutrino experiments 
(see also comments in [38]) unless the local density distribution is well known [39]. In view of 
this we will consider general (not spherically symmetric) density profile. 

Suppose neutrinos cross k layers of matter. The density jumps occur in the points x = Xi 
(i = !,-••, k — 1); x — xq and x = Xk are the points where neutrinos enter and leave the matter 
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Figure 5: The regeneration factor for the SNO charged current events integrated over different inter- 
vals of the observed kinetic energy as function of the zenith angle. We take Am 2 = 6.3 x 10~ 5 eV 2 , 
tan 2 9 = 0.4. 



correspondingly. Similarly to (60) we parameterize 9 m as 



Or. 



E sin 29 

Am 2 



k-l 



J2 ^Vi S(x - 



X i 



i=l 



where the jump of potential at the ith border between layers equals 

AVi = V( Xi + e) - V[x, - e), i — k-l. 



(6* 



(69) 



e is the infinitesimally small distance. Noting that the potential is zero for neutrinos before 
entering the Earth and after leaving the Earth, we define also 



AV = V{x ), AV k = -V{x k ). 
Plugging the potential jumps into (18) gives 



c(x fc ) 



Am 2 



(70) 



(71) 



i=i 



where 



dx——, i = 0, 



(72) 
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Figure 6: The same as in Fig. 5 for different intervals of averaging and the same middle kinetic energy 
T = 13 MeV. 

is the phase acquired from a given border i to the final point of the trajectory (detector). 

Now it is straightforward to compute the regeneration factor, and in the leading order in 
EV/Am 2 we obtain 



reg 



(u e \U(9 m (x k ))S(x k ,x )U\9 m (x ))U(9)\u 2 ) 



sin 2 9 



E sin 2 19 -^-^ 

— - — — > AV { COS 20;. 

Am 2 f^ Q 



(73) 



In the Appendix B, a direct computation of f reg for this case is given. Its result coincides 
with (73) in the leading order in EV/Am 2 . Using this formula one can easily reproduce (63) 
assuming a symmetric density profile and taking into account that <3?j = (0j— 0/c_i) for k = 2n—l 
and i < n. 

Using (73) it is easy to study averaging effects following the discussion in section 4.4. Struc- 
tures situated far from the detector have 0; ^> 1. So that, after averaging, remote structures 
do not produce significant effect. However, if the energy resolution is improved, effects of these 
remote structures can be large. This agrees with general consideration in [37]. 

Thus, we arrive at the following conclusion. If cos9z is large, small scale structures near 
the entering point can be taken effectively as a single layer. Furthermore, averaging over the 
energy makes contributions of these small structures to be unimportant. 

If cos 9z is small, we can not consider small scale structures near the entering point as a 
single layer. However uncertainties produced by these structures can be significantly reduced 
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if averaging is performed over broad energy interval, i.e. AE/E ~ 1. After averaging, the 
regeneration factor still shows an oscillatory behavior in the region of small cos 8z, and this 
effect is produced by contributions of the shells close to the detector [37]. 

5 Conclusion 

We have performed detailed analytic study of the LMA MSW conversion of the solar neutrinos. 
Our main result is the precise analytic formula for the survival probability which includes non- 
adiabatic corrections, averaging over the neutrino production region and the Earth regeneration 
effect. For the K component of the solar neutrino spectrum (K = pp,pep, Be, N, O, E, B, hep) 
it can be written as 

P K = \ + \{1- 5 k ) cos26 m (V K ) cos 26 - (1 - 5 K ) cos26 rn {V K )f Teg . (74) 

Here the correction due to averaging effect, 5k, is given in Eq. (35); the average values of 
matter potential in the production regions of K components, Vk, are defined in (30) and their 
numerical values are presented in the Table 1. The regeneration factor f reg is given in (62) for 
the symmetric density profile and in (73) for general asymmetric density profile. 

Effect of averaging over the neutrino production region in the Sun is reduced to specific 
value of the initial mixing angle in matter which should be taken for the average value of the 
potential, (P m = m (VK), and to the appearance of the correction 5k- We have compared the 
analytic results with the results of numerical computation and found that maximal deviation 
~ 1.8% happens for the hep neutrinos. For the boron neutrinos the precision is better than 
0.2%. 

We have obtained precise analytic formula for the regeneration effect in the Earth using the 
realistic density profile. We present simple derivation of this formula which uses the adiabatic 
perturbation theory. Performing also explicit calculations of the evolution in sequent layers 
we show that this derivation is correct. The analytic formula reproduces results of numerical 
computations with accuracy determined by r\ ~ 1 — 2%. 

Essentially the regeneration effect is the sum of contributions from different shells which are 
determined by jumps of the potential at the borders and by the adiabatic phase acquired inside 
the outer borders of the corresponding shells. The dependence of regeneration factor on the 
zenith angle can be understood in terms of interference of contributions from different borders. 

The derived analytical formula allows us to understand the effect of averaging over the neu- 
trino energy. Using the analytical formula we have considered effects of small scale structures 
(~ 10 km) of the Earth profile. These effects can be important for small values of cos#^. We 
stress that local "perturbations" of the density profile can produce sizable uncertainties in f reg . 
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Appendix A. Regeneration factor in a symmetric density 
profile 

Let us derive the expression for the regeneration factor by considering neutrino evolution in 
sequent layers of the Earth explicitly. First, we find the complete evolution matrix, S, in the 
basis of the mass eigenstates \I/ T = u 2 ): 

*(x f )=S(xf,x )*(x ). (75) 

As discussed in section 4.2, the adiabaticity violation effect within layers is suppressed by 
2E/(Am 2 R E ) ~ 1 — 2% in comparison with the leading order Earth matter effect (~ rj). 
Therefore, we neglect the adiabaticity violation within layers. 

1). In the case of neutrino propagation in one shell (one layer) we can simply project the 
adiabatic evolution matrix (11) obtained for the matter eigenstates on to the basis of the mass 
states. In the leading order in EV/Am 2 we find 



= u\e)u(e mR )s ad ~) u\e mR )u{6) = u(Ae m0 )s ad (%)u^(Ae m0 ) 



V2 2 J y ' VcosA^ m0 -smA0. 

e**° \ , EAV^nW ^ /0 1 



mO 
mO 



o °" e "Hi oj- (76) 

Here A9 rn0 = — # is the jump of the mixing angle at the surface of the Earth. We have 
used Eq. (50), and $ is defined in Eq. (46). 

2). In the case of two shells crossing, the neutrino encounters three layers (the outer shell is 
crossed twice). The evolution matrix can be similarly obtained by using the adiabatic evolution 
matrix (11) in each layer and by rotation from the matter eigenstates basis in the layer before 
the border to the basis after the border. As a result, we find 

S = U(A0,,^ (§, h) 4)5- (4, -§) ^(AM, 

where Si is the evolution matrix in the inner shell which has a form similar to Eq. (76) and it 
can be written as 

s, =u(Ae ml )s" d (^,-^u<(Ae ml ) 

Here A# m i is the jump of mixing angle on the border between the first and the second shells, 
and $i is defined in Eq. (46). 
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Combining the last two formulas we find to the order EV/ Am 2 

' L L" 



S =So 



2' 2 



= S 



ad 



2'-2) + £ 



K EAV sin 20 



i=0 



Am 2 



e -i*i _ e i*i 



1 

1 



3). Suppose the evolution matrix for j shells crossing (2j — 1 layers) equals 



S 



ad 



2' 2 
L L 



i=0 



1 

1 



Am 2 

Consider now the trajectory with j + 1 shells crossings. The evolution matrix is 
S = Sj+i 



(78) 



(79) 



= n 



U(A9 mi )S 



ad 



Li L. 



+i 



2 ' 2 



2 ' 



o 

n 

i=j-l 



s 



ad 



J i+l 



U\A0 mi ) 



,(80) 



where Ad m i defined in (48), is the jump of the mixing angle in matter at the border Ri. Sj is 
the evolution matrix in the central shell which can be written similarly to Eq. (77) as 



gad ( ^ _r£ ) + 



Am 2 



sin2#(e 



1 

1 



(81) 



$j is given in Eq. (46). After insertion into (80), the first term of (81) leads to Sj. The second 
term in (81) is already of the order EV/ Am 2 . Note that A9 m i is small, as is shown in (50). So, 
we can approximate U(A6 m i) by the unit matrix when plugging the second term in (81) into 
(80). As a result, we find 

1' 

1 



S j + 



£A^sin20 (e _^_ e , 



Am 2 



(82) 



Using then expression (79) for Sj, the formula (79) is immediately extended to the case of 
crossing j + 1 shells, thus accomplishing the proof. The result for the case of n shells crossing 
is 

or 

1 0, 



s = s 



ad 



L _L\ n ^ EAV iS m26 _ m 

Z Z/ i=0 



Am 2 



(83) 



Using (83) we obtain the expression for the regeneration factor in the leading order in 
EV/ Am 2 as 



reg 



■ a -<*o , Q n ^EAV iS m2e, 
sin 9e + cos 2^ 



2E sin 2 20 
Am 2 



i=0 
n-1 



Am 2 



- e" 



- sin 2 



sin $ AV i sin 



(84) 



i=0 



where $j is given in (46) and AVi is defined in (47). This expression coincides with (62) or (63) 
which have been obtained in section 4.3 using the adiabatic perturbation theory. 
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Appendix B. Regeneration in asymmetric density profile 

As in the section 4.5, we define 

A0» = 9 m (xi + e) - 9 m (xi - e), i = 0, • • • , k, (85) 

where A9 = 9 m (x ) — 9 and A9 k = 9 — 9 m (x k ), and the points of density jumps. 

xo and x k are the initial and final points of neutrino trajectory in matter. We will use the 
following expression: 

EAVj 

sin A9, t m — — ^ sin 29, cos A9 { « 1, (86) 
Am 2 

which is a good approximation in the leading order in EAV/ Am 2 . AV-i is given in (69). 
Neglecting the adiabaticity violation within each layer, we obtain the evolution matrix S as 

S= (l[U^A9 i )S ad (x i ,x^ 1 ))u\A9o). (87) 

S=fc ' 

S ad (x,xo) = S ad ((f)(x)) is the adiabatic evolution matrix given in (11). 
We approximate W(A9i) as 

£/+(A0O = (J J)+Q(Aft), (88) 



where 



Straightforward computation gives the following result in the leading order in EV/Am 2 
S = S ad (x k , x ) + Q(A9 k )S ad (x k , x ) + S ad (x k , x )Q(A9 ) 

+ J2S ad (x k ,x l )Q(A9 l )S ad (x l ,x ) 
i=i 

{ - sin A9-e i( - 2<i>i ~ <i> °} \ 

0j is defined in (72). Then, using (86), the regeneration factor can be directly computed in the 
first order in EV/Am 2 as 

k 2 

2 , 



f reg = |sin0e-^° -cos9J2 sin A^e i(2< ^ l,) 
E sin 2 29 



i=0 
,2 ofl 



sin'fl 



Am 



2 



^A\/,cos20 i . (91) 



i=0 



It coincides with (73). 
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Note added 



This note has been added on request of the referee. 

1) . After the present paper had appeared in the hep-ph archive [hep-ph/0404042], the 
preprint by Akhmedov et. al, [hep-ph/0404083], has been published in which the analytic 
integral formula is given for the the regeneration effect in the Earth in the three neutrino 
framework. In the first version of [hep-ph/0404083] the correct oscillation phase in this integral 
formula has been introduced on the "heuristic" basis: it does not follow from their perturbation 
theory. Correct integral formula (with the correct phase) has been derived for the first time in 
the paper by Ioannisian and Smirnov, [hep-ph/0404060]. Later in the Journal version JHEP 
0405 (2004) 057, Akhmedov et. al, have also presented derivation of correct phase. 

Let us now compare the results of papers [hep-ph/0404060], [hep-ph/0404083] with the 
results of present paper. (We will use the Akhmedov's et. al. results in the limit of zero 1-3 
mixing.) 

In Ioannisian and Smirnov paper hep-ph/0404060 and Akhmedov et. al. paper JHEP 0405 
(2004) 057 the integral formula has been obtained using the improved perturbation theory in 
the small parameter rj = 2EV(x)/Am 2 . In the present paper we use the adiabatic perturbation 
theory. It can be shown that in the lowest order in 77 both approaches coincide. Indeed, inserting 
expression for 6 m from (7) into (18) and performing integration by parts in Eq. (18) of the 
present paper one can derive the integral formula. 

This formula is valid for arbitrary density profile provided that the condition 77 1 is 
satisfied. Inserting the potential of PREM model (Eqs. (47, 53)) into the integral formula one 
can reproduce the result (63). However technically the use of formula (18) of the present paper 
is more convenient for the derivation of (63) since the derivative d6 m /dx gives 5- functions at 
the borders of layers and further integration becomes trivial. 

2) . The KamLAND collaboration has published recently results of improved measurements 
of oscillations on the basis of 766.3 ton-year exposure [T. Araki et. al, hep-ex/0406035]. In 
assumption of the CPT conservation, the global analysis of the solar neutrino data and Kam- 
LAND gives slightly (~ 10%) higher best fit value Am 2 = 8 • 10~ 5 eV 2 than it was before. The 
increase of Am 2 leads to the corresponding small decrease of the adiabaticity parameter 7 and 
the expansion parameter 77 for a given energy of neutrinos. Therefore increase in Am 2 (i) fur- 
ther improves the adiabatic perturbation theory and implies that the non-adiabatic correction 
for probability in the Sun is smaller; (ii) diminishes the Earth matter regeneration effect; (iii) 
makes our analytic formula for f reg preciser. Notice that the analytic study of this paper has a 
general character and does not rely on particular values of Am 2 . We use specific value of Am 2 
for illustration only. 
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